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This work presents a theoretical analysis on savings of energy during an endothermic reaction under 
microwave heating compared to conventional heating and shows the use of metal coated bounding 
surface to enhance the energy savings in otherwise low saving zones. Main thrust of this work is the 
quantification of energy savings for various probable microwave heating scenarios that may arise either 
due to varying reactor dimension (2L) over thin, intermediate and thick regimes or due to varying 
dielectric properties of the reactor. The analysis considers detailed transport equations in conjunction 
with Helmholtz equation for microwave propagation within a semiinfinite batch reactor. Simulations 
show that use of microwave can significantly save energy (as high as 60%) depending on reactor 
configuration. Simulations also show efficient use of metal coated bounding surface to enhance energy 
savings for reactors with 2L/A e ff = 0.5n-0.25, where n = 1, 2, 3... and / e ff is wavelength of microwave 
within the reactor. The enhancement is found to be 2 and 1.5 times at 2L/A e ff = 0.25 and 0.75, respec¬ 
tively. Various regions of efficient use of metal coated bounding surface for different microwave heating 
scenarios have been identified in a series of master curves. 

© 2013 Published by Elsevier Ltd. 


1. Introduction 

Energy consumption has been one of the important issues over 
past few decades triggering a series of research on finding faster 
and lower energy consuming technologies. Among several other 
possibilities like ultrasonic [1-16], electro-osmotic [17-20] and 
pulsed corona processing [21-26], microwave processing of ma¬ 
terials came out as an alternate viable process due to its faster 
processing and quick start up capabilities. Microwaves within fre¬ 
quency range of 300 MHz to 300 GHz interact with dipoles and can 
transfer energy volumetrically at molecular level in dielectric ma¬ 
terials. Volumetric supply of energy accelerates the progress of a 
process due to lower energy transfer time compared to cases where 
energy is supplied from surfaces. Over the years, potentials of mi¬ 
crowave processing in terms of process intensification and lower 
energy consumption have been proven in many different fields 
[27—39] with a major potential in the field of material synthesis, 
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where chemical reactions are reported to undergo dramatic ac¬ 
celeration under microwave radiations [40-60], 

Since reaction rates are strong functions of local temperature, a 
detailed knowledge about temperature profile inside the reacting 
domain is necessary to ascertain and quantify enhancement of re¬ 
action under microwave radiation. The situation becomes more 
demanding for endothermic reaction, where occurrence of reaction 
depends on local availability of heat and hence on heating pattern. 
It may be noted that microwave induced heating pattern may be 
highly complex with multiple local hot-spots within the domain 
depending on dimension and dielectric properties of the medium 
[61-65]. However, due to the difficulties associated with detailed 
measurement of electric field within the domain, microwave 
induced heating pattern has often been approximated by ad hoc 
models in literature. Use of such ad hoc models may hide some of 
the physics underlying the observed enhancement leading to 
speculation about athermal effect of microwave radiation on 
chemical reaction [66,67], Correspondingly, a few theoretical 
studies have been initiated to relate the enhanced reaction rate due 
to microwave heating effects [68-73], But, none of the earlier re¬ 
searchers considered detailed description of microwave propaga¬ 
tion via Helmholtz equation and thus could not quantify/analyze 
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Nomenclature 

T 

temperature, K 



To 

initial temperature, K 

Ca and ca dimensional (mol/m 3 ) and dimensionless reactant 

Z and z 

dimensional (m) and dimensionless spatial coordinate 


concentration 



Cao 

initial reactant concentration (mol/m 3 ) 

Greek symbols 

c a 

dimensionless average reactant concentration 

a e ff 

effective thermal conductivity (fe e fr/(pCp)eff) 

(. PCp)eff 

effective heat capacity (J/m 3 I<) 

P 

conduction number (k e ffTo/2I/o) 

Deff 

effective reactant diffusivity within packed bed (m 2 /s) 

AHr 

heat of reaction (J/mol) 

Dp eg 

effective penetration depth (m) 

Yr 

dimensionless activation energy (E/RT 0 ) 

E 

activation energy (kj/kmol) 

eo 

free space permittivity (Farad/m) 

Em 

electric field (V/m) 

k' eff 

effective dielectric constant of packed bed 

f 

frequency (Hz) 

k eff 

effective dielectric loss of packed bed 

fp 

l eff l2%D PeS 

■W 

effective wavelength of microwave within packed 

fw 

-WMo 


column, m 

Io 

intensity of incident radiation (W/m 2 ) 

^0 

wavelength within free space, m 

keff 

effective thermal conductivity of packed bed (W/m I<) 

X 

dimensionless time 

L 

half-width of packed column, m 

r D 

diffusion number {(j)a e ff/D e(f ) 

Np 

penetration number (2L/D Peff ) 

t /mw 

dimensionless reaction time under microwave heating 

Nr 

heat-reaction number (2IAHr 0R(Q o , T 0 )// 0 ) 

T /co„, 

dimensionless reaction time under conventional 

Nw 

wave number (2L//l e ff) 


heating 

QMwand q MW dimensional (W/m 3 ) and dimensionless absorbed 

e 

dimensionless temperature (k e ff(T— T 0 )/2t/o) 


power distribution 

<t> 

porosity 

<?MW 

dimensionless average absorbed power 

<S>R 

bulk Thiele modulus (4L 2 </>R(Qo, T 0 )/D eff C A0 ) 


time, s 




the entire spectrum of microwave heating possibilities and their 
influence on overall progress of reaction. 

Main objective of the current study is to fundamentally 
investigate the enhancement of reaction under various probable 
microwave heating scenarios (that can be encountered due to 
varying reactor dimension or dielectric properties of the reacting 
medium) and to show how the enhancement can be further 
improved by inserting a metallic coated surface at the other side 
of microwave incidence for certain specific cases. The reaction 
considered is a first order endothermic gas phase reaction 
occurring in a semiinfinite batch reactor in presence of susceptors. 
It may be noted that most of the common gases do not absorb 
microwave energy to significant extent and thus require suscep¬ 
tors to continue the reaction. Susceptors absorb microwave en¬ 
ergy and transfer them to the gaseous reacting medium. Reactions 
can also be carried out without the help of susceptors in cases of 
microwave absorbing reacting medium e.g. dielectric liquid re¬ 
actants. However, microwave propagation depends on reactant 
concentration within such system and change during the course 
of reaction depending on the specific material under consider¬ 
ation. Correspondingly, simulations of such microwave absorbing 
reacting system remain restricted to a specific material or a group 
of materials representing only a few specific microwave heating 
scenarios. This directed the current work to consider microwave 
non-absorbing gaseous reaction, which unleashes the limitations 
of material specific simulations and allows to perform a material 
invariant analysis considering all probable microwave propaga¬ 
tion scenarios. Nevertheless, the results presented here can be 
extended for qualitative estimates of microwave absorbing liquid 
phase reactions. 

The analysis is carried out via simulating progress of reaction 
from detailed mass and energy balance equations in conjunction 
with Helmholtz equation for microwave propagation. Simulations 
are carried out for microwave propagation with or without metal 
coated bounding surface in order to investigate the changes in 
reactor dynamics and probable enhancement of energy savings. 
Simulations are also carried out for reactions under equivalent 
conventional heating to quantify enhancement of reaction rate and 


corresponding reduction of energy consumption under microwave 
radiation in various probable scenarios of varying reactor di¬ 
mensions, presence of metallic coated bounding surface, changing 
susceptor properties, varying reaction rate etc. 

2. Mathematical modeling 

The reacting system under consideration is a semiinfinite 
batch reactor as shown in Fig. 1, where the entrapped gas un¬ 
dergoes a first order endothermic conversion given by 
A(g) -+ B(g). This work considers two possible configurations of 
heating the reactor via microwave as shown in Fig. 1(a) and (b). 
The right wall of the reactor is coated with a reflective (metallic) 
material in subplot (a), while the reactor does not contain any 
reflective coating in subplot (b). The reflective material 
completely reflects back the microwave and hence alters the 
resulting heating pattern in the former compared to the later. In 
both the cases, the reactor is exposed to microwave of intensity Io 
from the left end and packed with highly microwave absorbing 
materials or susceptors (e.g. SiC). Note that susceptors are 
required since gaseous species do not absorb microwave in 
general and thus cannot be heated under microwave radiation by 
themselves. The two configurations with and without reflective 
coating will be referred as “type 1” and “type 2" reactors, 
respectively. To compare the efficiencies of the two microwave 
heating strategies, their equivalent conventional heating with 
same heating rate are also considered as shown in Fig. 1(c). It 
may be noted that intensities of conventional heating corre¬ 
sponding to type 1 and type 2 configurations may differ from 
each other due to changes in microwave propagation. 

It is assumed that the reaction occurs only in the gas phase and 
the packing materials do not absorb any species. It is also assumed 
that material properties remain constant throughout the reaction, 
packings are homogenous and reactor walls are impermeable and 
thermally insulated. With further assumption of local thermal 
equilibrium and negligible convection in gas phase, evolution of 
reactant concentration (ca) and temperature (T) from the uniform 
initial state of Ca.o and To, respectively can be simulated from 
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following mass and energy balance equations and the associated 
boundary conditions: 

^ = D eS ^-<PR(C A ,T) (la) 


respectively, where subscripts g and p represent gas phase and 
packing material, respectively and D A and i p are molecular diffusivity 
of species A and tortuosity of packing, respectively. In addition, AHr 
is the heat of reaction and R[C A ,T) is the first order reaction rate 
satisfying Arrhenius temperature dependency given by 


(pC P ) en f t -kJ^ 2 + ^H R R(C A .T) 


f Qmw Microwave 
(0 Conventional 
(lb) 


Microwave 

L 

J Qm\n&Z Conventional 
- 1 

(lc) 

oz - L - § = °-S-° (ID 

These equations have to be solved in conjunction with Helmholtz 
equation for electric field (E m ) along with radiation boundary 
conditions for scattering of waves at reactor boundaries given in 
Appendix A (see Eqs. (A.1)—(A.6) for evaluation of absorbed power 
distribution within type 1 and type 2 reactors (Qmw) given by 

—SfcdXt)' 

where Ej), is the complex conjugate of E m , Ao and / e ff are wavelength 
of microwave with free space and reaction domain, respectively, D Peff 
is the penetration depth of microwave within reaction domain and 
E 0 is the incident field (see Appendix A). Other variables used in Eqs. 
(la)—(Id) are 21 as the reactor width, 0 as porosity of packing, 
Heff = D A <j>/Tp,(pC p ) eff = 0(/>C p ) g + ( 1 - 0)(pC p ) p andk e ff = 0k g + 
(1 - <j))kp as effective diffusivity of reactant, effective heat capacity 
and effective thermal conductivity of the packed reaction domain, 



R(Ca,T) = k R C A exp[-A] 


(3) 


where k R is the first order rate constant, E is the activation energy 
and R is the universal gas constant. 

In terms of the following dimensionless variables 

_ _ Z + L , = “efff a feeff(T-T 0 ) 

2 L ’ 4L 2 ’ 2 U 0 ’ 


2EQmw 

Io 


E o 


R(C A ,T) 

R(Qo-Jo) 


( 4 ) 


where a eff = k e{f /(pC p ) ef{ is the effective thermal diffusivity of the 
packed reaction domain, dimensionless form of governing mass 
balance, energy balance and Helmholtz equation and their associ¬ 
ated boundary conditions (Eqs. (la)—(Id) and Eqs. (A.1), (A.4) and 
(A.5a), (A.5b)) can be rewritten as 

(5a) 


00 

0T 


—2+N R r(c A ,ff) 


Qmw Microwave 
0 Conventional ’ 


(5b) 


d ^‘ + 4n 2 Ni(l + if p ) 2 Em = 0, 


(5c) 
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—a = o,-- = 


f 0 Microwave 
l Qmw Conventional ’ 


J w / w dz + ' £m = 2i ex P[-' TCN w/w] 


below a critical level (say r CT i t or 0 CT j t in dimensionless scale). This is 
to avoid unrealistic temperature drop during the simulation due to 
(5d) absorption of heat by reaction especially in no/low heating zones 
within the reactor. Reaction is resumed back once local tempera¬ 
ture increases to T cr i t . This causes the reaction rate to undergo a step 
jump given by 


@z = 



0 

i d E m 

2tt N w f w dz 


Type 1 
Type 2 


(5e) 


supplemented with the following expressions for dimensionless 
reaction rate (r(c A ,ff)) and dimensionless absorbed power (qMw) 

r(c A ,d) = c^exp ■ Qmw = E m E m *, (6) 


where dimensionless heat flux for conventional heating given by 
Qivrw — J z 1 = o < ?Mw(z)dz is the dimensionless average absorbed po¬ 
wer under microwave radiation. In this formulation, dimensionless 
parameters N w (wave number), f p and f w characterize wave prop¬ 
agation within the reactor and are defined as 


N w 


-W 

2n:Dp elf 


,/w 


7 e ff 

Ifc * 


(7) 


r(c A ,d) 


CA 6XP [ yR IT^\ for 19 ^ ^crit 
0 otherwise 


O) 


which causes the tracing of its dynamics to be numerically chal¬ 
lenging due to introduction of Dirac delta function in the Jacobian 
matrix [75], However, following the suggestion of an earlier work 
on tracing the phase change of pure materials [75], the numerical 
difficulties are efficiently circumvented by smoothing the step 
jump with a linear function within a superficial reaction zone of 
width 2AT around Tent as shown in Fig. 2. The modified rate 
expression with inclusion of superficial reaction zone can be writ¬ 
ten in terms of dimensionless variables as 


! 0 

rfa^crit + e) + 

r(c A ,8) 


for 6 < 0 crit - e 

for e crit -e< d < d crjt + e, 
for 6 > dan + e 

(10) 


where N w captures the effect of reactor thickness on wave propa¬ 
gation and dimensionless parameters f w and/ p capture the effect of 
dielectric properties on microwave induced heating. [It may be 
noted that the combination 2izN w f p — 2L/D Peg is the measure of 
reactor thickness compared to the penetration depth of microwave 
and is called penetration number or N p in literature [61,62,64] ]. The 
other four dimensionless numbers given by <E>r (bulk Thiele 
modulus), td (diffusion number or Lewis number), Nr (heat-reac¬ 
tion number) and [3 (conduction number) quantify relative mag¬ 
nitudes of various reaction and transport time scales within the 
reactor and are defined as 

= = (8) 

tR t a t a t AHg 

where t D - 4L 2 /D eff , t« = 4L 2 /a eff , t R = C A J4>R(C Aa ,T 0 ), 
tiH, = ((pC p ) eff T o )/(0AH R R(Q o ,To)) and t H = (2L(pC p ) eS T 0 )/I 0 
represent mass diffusion time, thermal diffusion time, reaction 
time, cooling time due to heat absorption by reaction and heating 
time, respectively. It may be noted that <h R and ((capture mass and 
heat transport resistances within the reactor, respectively, Nr 
measures relative supply of heat compared that required by reac¬ 
tion and id quantifies relative speed of mass and thermal diffusion 
within the reactor. 


3. Solution technique 

The governing equations (Eq. (5)) are spatially discretized via 
weak formulation of Galerkin finite element method using 50 
quadratic elements and time integrated using Crank-Nicholson 
scheme with a time step of At = 10 -4 . The resulting nonlinear 
algebraic equations are solved at each time step by Newton- 
Raphson method following the same procedure as laid out in an 
earlier article [74], It is worthwhile to mention here that simulation 
of endothermic reaction requires to stop the reaction locally and 
temporarily wherever and whenever local temperature reduces 


where e - k e ffAT/2L/ 0 is half width of the superficial reaction zone in 
dimensionless temperature scale. Convergence of numerical sim¬ 
ulations with respect to the width of superficial reaction zone are 
shown in Fig. 3, which shows variations of overall reactant deple¬ 
tion rate (1 - c A ) as r is varied from 10 -3 to 10 -6 (dotted, dash and 
solid lines and circles). Simulations in Fig. 3 are carried out for low 
intensity conventional heating corresponding to type 1 and type 2 
configurations at N w — 0.06, f w = 0.1, f p — 0.1 (q MW = 0.002 and 
0.141, respectively). It may be noted that unrealistic temperature 
drop during the simulation become more probable for surface 
heating (at regions far away from heat source) than volumetric 
heating of microwave especially for low intensity of q MW . As a 
result, simulations require much frequent use of superficial reac¬ 
tion zone (Eq. (10)) in cases of low intensity of conventional heat¬ 
ing. Correspondingly, convergence of numerical simulations in the 
two severe cases presented in Fig. 3 guarantees convergence for all 
other cases. Other parameters in Fig. 3 correspond to <1> R = 10, 
t d — 1, N r — 0.1, /? = 0.1 and jr — 10. Fig. 3 shows that numerical 
simulations become invariant of the width of superficial reaction 


r(c a ,t) 
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<1mw= 0 - 002 <J yP e !) 



10' 2 10' 1 1 

Dimensionless time (x) 


9w= 0 - 141 (Type 2) 



lCT 2 10 4 

Dimensionless time (x) 


Fig. 3. Convergence of numerical simulation with respect to the width of superficial reaction zone (c) shown in terms of 1 - c„vs t at e = 10 -3 (dotted lines), 10 -4 (dashed lines), 
10 5 (solid lines) and 10“ 6 (circles) for two severe cases of slow surface heating corresponding to type 1 and type 2 configurations at N w = 0.06 ,f p = 0.1 and/ w = 0.1 with q MW = 
0.002 and 0.141, respectively with other parameters as t d = 1, N R = 0.1, 0 = 0.1, <t> R = 10 and y R = 10. 


zone at e = 10 -5 even for these severe situations. Hence, we will 
use e — 10 -5 for all the simulations. It may be noted that we have 
assumed r cr it = To in Fig. 3 and will use the same in all forthcoming 
illustrations. 

Starting from ca = 1 and 6 = 0 at x = 0, all simulations are 
stopped when conversion reaches 90% across the entire reactor 
corresponding to max(c/i) = 0.1. Time required for this conversion 
will be referred as reaction time and denoted by iy MW and ij Qinv for 
microwave and conventional heating, respectively. Simulations are 
carried out for different reactor dimensions over the entire range of 
N w to replicate various possible microwave heating scenarios and 
their alteration due to metallic coating. For all the cases, reactor 
dynamics are simulated and compared for two limiting cases of 
weak and strong mass transfer limitations within the reactor via 
<I>r = 0.1 and 10, respectively with fixed relative heat supply as well 
as fixed heat transfer limitations at Nr = 0.1 and ft = 0.1, respec¬ 
tively. ft may be noted that Nr<^\ ensures enough heat supply to 
the reactor compared to the heat requirement for reaction. On the 
other hand, higher thermal diffusion time compared to heating 
time at /?«! allows to simulate distinct dynamics of microwave and 
conventional heating [74], It is important to note that heating of 
reactor becomes invariant of the supply pattern for /?»1, where 
supplied heat spreads uniformly across the reactor instantaneously. 
Other parameters are selected to be t d = 1 (to avoid preferential 
diffusion of either heat of mass) and jr - 10. 


4. Model validation 

In absence of any appropriate experimental or numerical data in 
literature, we validate the present simulations by comparing them 
with the analytical solutions which can be obtained in the limit of 
jr —> 0. In this limit of jr —► 0, the nonlinear mass and energy 
balance equations (Eqs. (5a) and (5b)) reduce to the following linear 
ones: 




(lla) 


6 e__^e f q mw Microwave 

9x a? + R A ~ \ 0 Conventional ’ 


(lib) 


along with the boundary conditions given by Eqs. (5d) and (5e) and 
the initial conditions given by c A = 1 and 6 = 0 at x = 0. Using Finite 
Fourier Transformations (FFT), closed form solutions of the above 
equations are obtained as 


c ^* e *p [-$']’ 


(12a) 


0(z,t)=xq MW -fe( l-exp[-^x]) 

“ 1 - exp[->. n x] , s f (<?MW,Wn) Microwave 

+ E--P- -w n (z)< . 

n=i A n L v 2q MW Conventional, 

(12b) 

where A n = nV is the nth eigenvalue, w n (z) is the nth eigenfunc¬ 
tion given by 

"'■«-{ Eos ,«z S>1°. < 12c > 

and (q MW , w„) is the inner product defined as 

(< 7 mw,w„) = J q MW (z)w„(z)dz. (12d) 


The expressions for (qmw, w n) and q MW are given in Appendix B 
(Eqs. (B.4) and (B.5)). 

To compare numerical results, reactor dynamics in the limit of 
jr —> 0 are simulated from Eq. (5) by using jr = 0.001 with other 
parameters as 4>r = 10, Nr — 0.1, td = 1.0 and ft = 0.1 (same as in 
Fig. 3). Resulting concentration and temperature profiles are 
compared with those obtained from the analytical expressions 
(Eqs. (12a) and (12b), respectively) for three different reactor di¬ 
mensions of N w — 0.25, 0.5 and 0.75 with f p — 0.1 and f w = 0.1 as 
shown in Fig. 4(a)-(b). These three reactor dimensions correspond 
to various probable stiffness of microwave heating patterns (as 
reflected in the corresponding temperature profiles) and also 
correspond to a wide range of conventional heating intensities 
(q M w = 0.949, 0.114 and 0.6 for type 1 configuration and 
q MW = 0.059,0.504 and 0.158, for type 2 configuration at N w = 0.25, 
0.5 and 0.75, respectively). Fig. 4 shows the comparisons for both 
microwave (MW) and conventional (Conv.) heatings, where bullets 
represent numerical simulations and continuous lines represent 
analytical solutions with darker and lighter shades indicating type 
1 and type 2 configurations, respectively. The spatial profiles in 
Fig. 4 are shown at x = 0.15, which is half the time required for 
completion of 90% conversion [Note: analytical solution predicts 
required reaction time as x = 0.2303, while simulations predict 
required reaction time to be x = 0.2307], It may be noted that in the 
limit of jr —> 0, reactant concentration and its consumption be¬ 
comes invariant of heating mode and hence invariant of spatial 
location within the reactor (as reflected in Eq. (12a)). Accordingly, 
numerical simulations predict uniform concentration for both mi¬ 
crowave and conventional heating irrespective of reactor dimen¬ 
sion as shown in Fig. 4(a). In contrast to concentration, temperature 
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(a) 



Dimensionless position, z 


Analytical: Continuous lines 

-Type 1 

-Type 2 

Numerical: Bullets 

• • • Type 1 

• « « Type 2 



Fig. 4. Comparison between simulated results (bullets) and analytical solutions (continuous lines) for type 1 and type 2 configurations (represented by darker and lighter shades, 
respectively) in case of three representative reactor dimensions of N w = 0.25,0.5 and 0.75 with/ p = 0.1, f w = 0.1, N R = 0.1, z D = 1, <Fr = 10, fl = 0.1 and y R = 0.001. Subplot (a) shows 
comparison of concentration profiles at t = 0.15 and subplot (b) shows comparison of temperature profiles for microwave heating (MW) and corresponding conventional heating 
(Conv.) at t = 0.15. 


profiles depend on heating patterns and thus vary with heating 
mode, reactor configuration and reactor dimension as shown in 
Fig. 4(b). Nevertheless, numerical simulations are in excellent 
agreement with analytical solutions in all the cases showing the 
accuracy of the solutions for different probable heating scenarios. 

5. Microwave power absorption and influence of metallic 
coating 

Influence of metallic coating on microwave power absorption 
are measured by two factors: (i) changes in the spatial distribution 
of absorbed power (qmw) and (ii) variations in the intensity of q MW 
or surface heat flux. Spatial distribution of absorbed power de¬ 
termines the homogeneity of microwave heating, which in 
conjunction with intensity of conventional heating (q MW ) play a 
critical role in determining the efficiency of microwave heating. 
Influence of metallic coating on the above two factors are shown in 
Fig. 5, where the main plot shows the variations of q MW and the 
representative insets show the variations of absorbed power dis¬ 
tributions from type 1 (solid lines) to type 2 (dotted lines) config¬ 
urations for various probable microwave heating scenarios 
spanning the range of reactor dimension given by 0.05 < N w < 10. 
The results in Fig. 5 are simulated using f p — 0.1 and f w — 0.1. 

The insets of Fig. 5 show that entire spectrum of power ab¬ 
sorption patterns can be categorized in three distinct classes of 
uniform (N w = 0.06), (ii) oscillatory ( N w — 0.5 and 0.75) and (iii) 
exponentially attenuated ( N p — 3) profiles. These three patterns 
appear in chronological order as reactor dimension increases from 
much below the wavelength of microwave (N w < 0.1) to much 
greater than penetration depth of microwave (N p > 3) and divide 
the entire range of reactor dimensions in thin ( N w < 0.1), interme¬ 
diate (0.1 < N w and N p < 3) and (iii) thick ( N p > 3) regimes, 
respectively [62], Thin regime refers to the entire range of reactor 
dimensions below N w = 0.1, where microwave power absorption 
occurs uniformly irrespective of metallic coating as seen in the inset 
for N w = 0.06. However, metallic coating significantly reduces mi¬ 
crowave power absorption in thin regime as can be seen from the 
magnitudes of q MW for N w < 0.1 in Fig. 5. It may be noted that 


microwave power absorption in type 1 configuration becomes 
insignificantly small as N w decreases beyond 0.05. 

Thick regime refers to the reactors with dimensions much 
greater than penetration depth of microwave ( N p >3), where mi¬ 
crowave losses most its energy before reaching the other end of the 
reactor. Correspondingly, power absorption in thick regime occurs 
only near the exposed surface and remain unaltered by the pres¬ 
ence of metallic coating as seen in Fig. 5 for N p > 3. It may be seen 
from the inset for N p = 3 that absorbed power attenuates expo¬ 
nentially from the exposed surface in thick regime. Power absorp¬ 
tion regime further shrinks with increasing reactor dimension and 
microwave heating approaches surface heating in the limit of 
Np -► OO. 

Intermediate regime refers to the reactors in between thin and 
thick regimes (0.1 < N w < 3/27t/ p ), where reactor dimensions are 
neither much lower than wavelength of microwave nor much 
greater than penetration depth of microwave. In this regime, mi¬ 
crowave power absorption is strongly dependent on the nature of 
scattering from the reactor walls. Scattering of microwaves form 
various traveling waves within the reactor and interferences be¬ 
tween them determine the resultant power absorption character¬ 
istics depending on phase-space of individual wave [62], Since 
metallic coating alters scattering by completely reflecting micro- 
waves from the right wall, resulting traveling waves in type 1 
configuration also differ significantly from those in type 2 config¬ 
uration. As a result, intermediate regime exhibits most prominent 
influence of metallic coating as seen in Fig. 5. 

Most important characteristic of intermediate regime is the 
oscillatory variations of q MW with reactor dimension, which result 
from constructive and destructive interferences between traveling 
waves. Constructive interference enhances power absorption, 
while destructive interference suppresses power absorption. In 
intermediate regime, constructive interferences occur at periodic 
intervals as reactor dimension increases from 0.1 to 3/2 tt/ p and 
every two consecutive constructive interferences are accompanied 
by a destructive interference in between [62], Hence, q MW un¬ 
dergoes periodic enhancement and suppression over the range of 
intermediate regime, which is reflected in Fig. 5 as a series of 
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10 1 1 

Penetration number ( N p ) 


Fig. 5. Influence of metallic coating on over the entire range of thin to thick regimes (0.05 <N W < 10) as well as influence of metallic coating on the power absorption profiles 
[<Jmw(z)1 at representative reactor dimensions (N„ or N p ), where solid and dotted lines represent reactors with (type 1) and without (type 2) metallic coating, respectively and 
parameters f w and f p are taken to be 0.1. 


maxima and minima of q MW for both type 1 and type 2 reactors. 
Enhancement of q MW at locations of its maxima are also called 
resonances. Fig. 5 shows that locations of resonances depend 
strongly on reactor configuration, where resonances for type 1 
configuration occur at the locations of minima for type 2 configu¬ 
ration and vice versa. Resonances occur at N w = 0.5n, n = 1,2... for 
type 2 configuration, where type 1 configuration leads to minima of 
q MW . Similarly, resonances for type 1 configuration occur at 
N w = 0.5n-0.25, n = 1,2..., where q MW exhibits minima for type 2 
configuration (see Fig. 5). Amplitudes of resonances are also much 
higher in type 1 configuration than type 2. 

Similar to q M w. power absorption patterns of intermediate 
regime are also strongly influenced by metallic coating, where type 
1 and type 2 configurations lead to completely opposite absorbed 
power distributions as illustrated in Fig. 5 for N w = 0.25, 0.5 and 
0.75. Typical feature of intermediate regime is the presence of local 
maxima in absorbed power distributions due to occurrence of 
spatial resonances, where larger reactor favors more spatial reso¬ 
nances as can be seen from the insets of Fig. 5 for N w = 0.25,0.5 and 
0.75. In all the cases, locations of spatial maxima in type 1 config¬ 
uration correspond to locations of local minima of absorbed power 
in type 2 configuration and vice versa. 

6. Alteration of reactor dynamics under microwave heating 
and influence of metallic coating 

Alteration of reactor dynamics under thin, intermediate and thick 
regimes of microwave heating for type 1 and type 2 configurations 
are illustrated in Figs. 6-10 for the representative reactor di¬ 
mensions corresponding to N w = 0.06 (Fig. 6), 0.25 (Fig. 7), 0.5 
(Fig. 8), 0.75 (Fig. 9) and N p = 3 (Fig. 10). In each case, influences of 
microwave heating are illustrated for two limiting cases of weak 
and strong mass transfer limitations at 4 >r = 0.1 and 10 in subplots 


(a) and (b), respectively. It may be noted that reactor dynamics 
under microwave heating are shown at the top panel in each 
subplot of Figs. 6—10 and are denoted by “MW”. Reactor dynamics 
under conventional heating corresponding to type 1 and type 2 
configurations of various N w 's and N p are shown at the bottom 
panel in each subplot of Figs. 6-10 and are denoted by “Conv”. In all 
the cases, parameters f p and f w are considered to be 0.1 and type 1 
and type 2 configurations are represented by solid and dashed 
lines, respectively. Here, reactor dynamics are compared in terms of 
spatial non-uniformities of reaction pattern ( N r t{ca,6 vs z) and 
reactant concentration {ca vs z ) within the reactor as well as overall 
reactant depletion rate (c^vs t). In each figure, spatial non¬ 
uniformities within the reactor are shown at two intermediate 
time levels given by t// 3 (darker shade) and 2iy/3 (lighter shade). 
Power absorption profiles along with magnitude of q MW for type 1 
and type 2 reactor configurations of Figs. 6-10 are also shown at 
top of each figure for ease of illustrations. 

Alteration of reaction dynamics under microwave heating is 
purely driven by individual power absorption patterns of type 1 
and type 2 configurations in thin, intermediate and thick regimes. 
Conventional or surface heating always results in localized reaction 
from the hot left surface, where intensity of reaction is purely 
determined by the surface heat flux or q MW . On the other hand, 
reaction occurs volumetrically under microwave heating but with 
varying spatial profiles depending on the regime of operation and 
reactor configuration. Reaction occurs uniformly under thin regime 
of microwave heating (N w < 0.1) as reflected in almost uniform 
spatial profiles of N R r(CA,ff ) at N w — 0.06 for both type 1 and type 2 
configurations. In contrast, reaction occurs only near the exposed 
left surface (similar to conventional heating) in thick regime of 
microwave heating as seen in the spatial profiles of N R r(CA,0) for 
N p = 3 in Fig. 10. In between these two regimes, reaction patterns 
under microwave heating vary with reactor dimension as well as 
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Fig. 6. Spatial dynamics ( N R t(c A ,8 ) vs z and c A vs z shown at iy/3 and 2t^3 in darker and lighter shades, respectively) and overall progress (c^vs t) of type 1 (solid lines) and type 2 
(dotted lines) reactors under microwave heating (MW) along with those under their equivalent convectional heating (Conv.) in thin regime corresponding to N w = 0.06, f w = 0.1 and 
f p = 0.1 with T D = 1.0, N r = 0.1, (3 = 0.1, y R = 10 and $ R = 0.1 and 10 in subplots (a) and (b), respectively. 
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Fig. 9. Spatial dynamics (N R r(c A ,0) vs z and c A vs z shown at ty/3 and 2q/3 in darker and lighter shades, respectively) and overall progress (c A vs t) of type 1 (solid lines) and type 2 
(dotted lines) reactors under microwave heating (MW) along with those under their equivalent convectional heating (Conv.) in intermediate regime corresponding to N w = 0.75, 
f w = 0.1 and f p = 0.1 with t d = 1.0, N R = 0.1, 0 = 0.1, y R = 10 and = 0.1 and 10 in subplots (a) and (b), respectively. 
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Fig. 10. Spatial dynamics (JV R r(cA,0) vs z and c A vs z shown at t// 3 and 2q/3 in darker and lighter shades, respectively) and overall progress (c A vs i) of type 1 (solid lines) and type 2 
(dotted lines) reactors under microwave heating (MW) along with those under their equivalent convectional heating (Conv.) in thick regime corresponding to N p = 3 ,f w = 0.1 and 
f p = 0.1 with t d = 1.0, N r = 0.1, 0 = 0.1, y R = 10 and $ K = 0.1 and 10 in subplots (a) and (b), respectively. 


configuration but remain bound within uniform reaction of thin 
regime and surface reaction of thick regime. Reactions mostly occur 
volumetrically under intermediate regime of microwave heating 
except near the boundaries of thin and thick regimes. It is inter¬ 
esting to note that similar to thick regime, reaction under micro- 
wave heating also becomes highly one sided within a small region 
near the boundary of thin regime around N w = 0.25, where reaction 
is driven from the reactor boundary as seen in Fig. 7 for N w = 0.25. 

Special characteristics of reaction under intermediate regime of 
microwave heating is presence of multiple reaction fronts at the 
locations of spatial resonances especially for larger reactors 
(N w > 0.5). Reactor with dimensions N w < 0.5 leads to a single 
reaction front irrespective of reactor configuration as seen in Fig. 7 
for N w = 0.25. However, location of reaction front changes from left 
to right surface as reactor configuration changes from type 1 to type 
2 at N w — 0.25 (see Fig. 7). Location of reaction front under mi¬ 
crowave heating also changes from the center to both the surfaces 
as reactor configuration is changed from type 1 to type 2 at 
N w — 0.5. This interestingly introduces two reaction fronts in type 2 
configuration at N w — 0.5 in contrast to the single reaction front in 
type 1 configuration (see Fig. 8). Reactors with N w > 0.75 lead to 
multiple reaction fronts irrespective of reactor configurations, but 
reaction fronts of type 1 configuration always appear in the loca¬ 
tions of minimum reaction zones of type 2 configuration as illus¬ 
trated in Fig. 9 for N w — 0.75. 

Common feature of all the reaction patterns is homogeneity of 
reaction under microwave heating compared to their conventional 
counterparts except for very thick reactors with N p -> °o. Corre¬ 
spondingly, both type 1 and type 2 reactors exhibit homogenized 
concentration profiles under microwave heating in presence of 
mass transfer limitations at <E>r = 10 as illustrated in subplot (b) of 


Figs. 6—10. Mass transfer limitations exist when diffusion of reac¬ 
tant is slower than its consumption rate or equivalently <&r > 1. In 
such cases, reactant builds up in no/low reaction zones of the 
reactor according to the non-uniformity of reaction pattern. Volu¬ 
metric power absorption significantly enhances the uniformity of 
reaction from conventional to microwave heating in thin and in¬ 
termediate regimes (except around N w ~ 0.25). As a result, reactors 
with N p < 3 and N w * 0.25 show prominent homogenization of 
concentration gradients under microwave heating compared to 
conventional heating as shown in Figs. 6(b), 8(b) and 9(b) for 
N w — 0.06,0.5 and 0.75, respectively. On the other hand, highly one 
sided reaction from the surface significantly reduces homogeniza¬ 
tion ability of microwave heating around N w ~ 0.25 as well as in 
thick regime as illustrated in Figs. (7b) and (10b) for N w = 0.25 and 
N p = 3, respectively. 

Homogenization ability of microwave heating also depends on 
the degree of localization in its conventional counterpart and thus 
on the intensity of q MW . Reactor dynamics remain homogenized 
(uniform) even under conventional heating if intensity of q MW re¬ 
mains sufficiently low such that resulting spatial gradients can be 
easily homogenized by diffusion. For such low intensity of q MW , 
reactor dynamics become invariant of heating pattern as is the case 
for type 1 configuration at N w - 0.06 (q MW = 0.0016), where mi¬ 
crowave and conventional heating lead to almost same reactor 
dynamics at = 0.1 as well as at 'I> K = 10 (see the profiles in solid 
lines). Influence of heating pattern increases with increasing q MW 
and becomes prominent at sufficiently high intensity of q MW , 
where spatial gradients build up in presence of mass transfer lim¬ 
itations (<E>r > 1) according to the non-uniformity of individual 
heating pattern. Correspondingly, entire intermediate regime of 
microwave heating with significantly high q MW as well as with 
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volumetric power absorption (except for the region around 
N w ~ 0.25) leads to much homogenized reactor dynamics 
compared to conventional heating at <f> R > 1 irrespective of reactor 
configuration as seen in Figs. (8b) and (9b) for N w = 0.5 and 0.75. 
Microwave heating also exhibits homogenization in thin regime for 
type 2 configuration, where q MW is not insignificantly small (see 
profiles of Fig. 6(b) in dashed lines). Highly one sided power ab¬ 
sorption near the region around N w = 0.25 obviously decreases the 
homogenization ability of microwave heating in spite of extremely 
high q MW for type 1 configuration and microwave heating also 
leads to significant concentration gradients within the reactor as 
seen in Fig. 7(b) for N w = 0.25. Similarly, microwave and conven¬ 
tional heating lead to almost similar spatial gradients at N p = 3 as 
shown in Fig. 10. 

Superior homogenization ability of microwave heating over 
conventional heating becomes irrelevant in the limit of ( I>r -c 1, 
where faster diffusion than consumption of reactant homogenizes 
the concentration gradients irrespective of heating patterns. 
However, magnitude of <I>r required to reach this limit depends on 
the non-uniformity of heating/reaction pattern and hence on q MW 
for conventional heating. Faster diffusion of reactant at 3 >r = 0.1 is 
observed to completely homogenize the concentration gradients 
for all the cases of microwave heating in volumetric power ab¬ 
sorption regimes (i.e. for N w = 0.06,0.5 and 0.75 (Figs. 6(a), 8(a) and 
9(a), respectively). Concentration gradients within conventional 
reactors are also completely homogenized at <I>r = 0.1 for low to 
moderately high q MW corresponding to (i) both type 1 and type 2 
configurations at N w = 0.06 [Fig. 6(a)], (ii) type 2 configuration at 
N w — 0.25 and 0.75 (Figs. 7(a) and 9(a), respectively) and (iii) type 1 
configuration at N w — 0.5 (Fig. 8(a)). However, concentration gra¬ 
dients within conventional reactors are not completely homoge¬ 
nized at fI>R = 0.1 for cases of significantly high q MW corresponding 
to (i) type 1 configuration at N w = 0.25 and 0.75 and (ii) type 2 
configuration at N w — 0.5. In such cases, conventional reactors 
exhibit slight concentration gradients at <f> R = 0.1, which disappear 
with further decrease of Thiele modulus. For the same reason, slight 
concentration gradients are also observed at 4> R = 0.1 for both 
microwave and conventional heating at N p = 3 (see Fig. 10(a)). 

7. Savings of energy under microwave heating: influence of 
metallic coating 

Since same heating rate is maintained, variations of energy 
consumption in microwave heating result from the differences in 
total processing time compared to conventional heating. For an 
endothermic reaction, changes in processing time are mainly 
caused by the diffusion time required by the reactant to reach the 
heating zone/zones from the colder regions of the reactor prior to 
the occurrence of reaction. It is obvious that volumetric supply of 
heat reduces diffusion time in microwave heating compared to the 
highly one sided conventional heating, where reactants from cold 
right end of the reactor have to diffuse the entire reactor width to 
avail the heat and react. As a result, microwave heating reduces the 
processing time as seen in all the c^vs t diagrams of Figs. 6—10 
especially for <P R = 10. Associated reduction of energy consump¬ 
tion can be quantified as 



where ij M w determines percentage savings of energy in microwave 
heating compared to conventional heating. Extent of energy sav¬ 
ings obviously depends on the homogenization ability of micro- 
wave heating and thus changes with reactor dimension (N w ) as well 
as reactor configuration as reflected in the c^vs x diagrams of 


Figs. 6—10. These changes are quantified in Fig. 11(a)— (b), which 
shows the variations of ?/mw with N w (solid lines) over the entire 
range of thin to thick regimes at various Thiele modulus of <I>r = 0.1, 
1 and 10 for type 1 [subplot (a)] and type 2 [subplot (b)[ configu¬ 
rations. Corresponding q MW vs N w diagrams are also shown by 
dotted lines in Fig. 11(a) and (b). As expected, savings of energy 
under microwave heating increases with increasing mass transfer 
limitations (Thiele modulus) and becomes prominent for <f> R > 1, 
where microwave heating can potentially reduce energy con¬ 
sumption as reflected in the magnitude of ?)mw- Fig-11 shows the 
role of reactor dimension (N w ) to maximize the utilization of mi¬ 
crowave energy for a given reactor configuration. Fig. 11 also shows 
the role of metallic coating to further enhance the energy savings at 
specific reactor dimension. It may be noted that savings of energy 
can be significantly enhanced at N w — 0.5n—0.25 by selecting re¬ 
actors with metallic coating as in type 1 configuration. On the other 
hand, metallic coating should be avoided at N w = 0.5n as well as in 
thin regime (N w < 0.1). Obviously, savings of energy remain unin¬ 
fluenced by metallic coating in thick regime. 

It follows from the discussion on Figs. 6—10 that two factors 
given by (i) variations of the homogeneity of power absorption 
pattern and (ii) changes in the magnitude of q MW influence the 
variations of ijmw as N w varies from thin to thick regimes. It may be 
noted from Fig. 5 that changes in the homogeneity of power ab¬ 
sorption with N w are much lower compared to the changes of q MW 
over the entire volumetric power absorption regime of N p < 3. 
Correspondingly, ijmw mostly follows the variations of q MW for 
N p < 3 in both Fig. 11(a) and (b) except for a very small region near 
N w — 0.5 in case of type 1 configuration as shown by shading in 
Fig. 11(a). It may be noted that as N w increases within the shaded 
region of Fig. 11(a), ijmw increases while q MW undergoes a reduc¬ 
tion. This anomaly interestingly leads to an additional peak in ijmw 
at N w = 0.52 for type 1 configuration in presence of mass transfer 
limitations as may be observed from Fig. 11(a) for <f> R = 1 and 10. 

Anomalous variations of q Mw within the shaded region of 
Fig. 11(a) are due to the changes in homogeneity of absorbed power, 



f p = 0.1, z D = 1.0, N R = 0.1, (3 = 0.1,= 10 and varying Thiele modulus as 4> R = 0.1,1 and 
10. Corresponding q MW vs N w diagrams are also shown in dotted lines (corresponding 
scale is shown at right axis). 
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Fig. 12. Influence of/ p on variation patterns of and ij M w with N w for type 1 (solid 
lines) and type 2 (dotted lines) configurations at f w = 0.1 and = 10. 


which undergoes a significant transformation from N w — 0.1 to 0.5 
irrespective of the reactor configuration as shown in the insets of 
Fig. 5. It may be noted that homogeneity of power absorption re¬ 
duces drastically from N w = 0.1 to N w = 0.25 and then reverts back 
as N w increases to 0.5. Increasing homogeneity of absorbed power 
causes t/mw to increase from N w = 0.4 to N w ~ 0.5 (more specif¬ 
ically from N w = 0.38 to 0.52 at 4> R = 10 and from N w = 0.45 to 0.52 
at <b R = 1) in spite of decreasing q MW for type 1 configuration. It 
may be noted that reduction of q MW are not significant within the 
shaded region. On the other hand, changes of q MW are exceptionally 
high within the range of 0.15<N W <0.4 for type 1 configuration 
and hence rj M w follows the same regardless of the changes of ho¬ 
mogeneity of power absorption. This causes t/mw to increase from 
N w = 0.1 to 0.25 for type 1 configuration in spite of decreasing 
homogeneity of absorbed power. Similarly, t/mw decreases from 
N w — 0.25 to N w ~ 0.4 in Fig. 11(a) although homogeneity of power 
absorption increases over that range. However, changes in the 
homogeneity of power absorption suppress the variations of t/ M w at 
N w = 0.25 compared to the corresponding resonating peak of q MW 
as may be observed from Fig. 11(a). Also, since mass transfer 


limitations bring out the effect of varying homogeneity of absorbed 
power distribution, the region of anomalous variations of t/mw also 
shrinks with decreasing <b R and disappear at 4> R = 0.1. In contrast to 
type 1 configuration, influence of the variations of q MW corroborate 
with influence of changing homogeneity of power absorption over 
the entire range of 0.1 < N w < 0.5 for type 2 configuration and ijmw 
shows one-to-one correspondence with q MW in Fig. 11(b) for the 
entire range of N p < 3 irrespective of Thiele modulus. 

8. Effect of f p on influence of metallic coating 

Influence of/ p on q MW over the entire range of thin to thick re¬ 
gimes (0.05 < N w < 10) and associated changes in tjmw at 4> R = 10 
for type 1 and type 2 configurations are shown in Fig. 12 by 
increasing/p from 0.1 to its upper bound of 1 [Note that f p is bound 
within 0—1 according to Eq. (A.3)]. In Fig. 12, solid and dotted lines 
represent type 1 and type 2 configurations, respectively. Here, mass 
transfer limited situations of <b R = 10 is selected to prominently 
bring out the effect of/ p on ijmw- All the simulations in Fig. 12 are 
performed at f w — 0.1 with all other parameters to be same as in 
previous figures. Fig. 12 shows that resonances of power absorption 
are strongly affected by/ p , where resonating peaks of q MW gradu¬ 
ally disappear as / p increases beyond 0.5 (or 1.5/rc more precisely 
[62]). It may be noted that most of the resonating peaks of q MW , 
except the first few, disappear even at f p — 0.3 for both the con¬ 
figurations. With suppression of resonances, intermediate regime 
also shrinks and causes its volumetric power absorption to trans¬ 
form to the highly one sided distribution of thick regime as f p in¬ 
creases from 0.1 to >1.5/tc as illustrated in Fig. 13. Fig. 13 shows the 
transformation of power absorption profiles in type 1 and type 2 
configurations as f p increases from 0.1 (solid line) to 0.5 (dashed 
line) and 1 (dotted line) at N w = 0.06, 0.25,0.5 and 0.75 for/ w = 0.1. 
As power absorption transforms towards the highly one sided 
distribution of thick regime, influence of metallic coating also 
gradually decreases with increasing^, as shown in Fig. 13. It may be 
noted that power absorption profiles of N w = 0.5 and 0.75 become 
invariant of the reactor configuration beyond f p = 0.5. Corre¬ 
spondingly, region of influence of reactor configuration on q MW 
shrinks considerably from N w <5 atf p = 0.1 to N w < 1.5,0.9 and 0.5 
at f p — 0.3,0.5 and 1, respectively (see Fig. 12). As a result, t/mw also 
shows a decreasing region of influence of reactor configuration 
from N w < 4 to N w < 1.5, 0.6 and 0.5 as f p increases from/ p = 0.1 to 
f p = 0.3, 0.5 and 1 respectively. 

Transformation to highly one sided power absorption profiles 
reduces the homogeneity of microwave heating as f p increases from 
0.1 to 1 for N w > 0.5, Correspondingly, Fig. 12 shows a significant 
reduction of rj M w from f p = 0.1 to 1 for both type 1 and type 2 
configurations over the entire range of N w > 0.5. Over this range of 



Fig. 13. Variations of power absorption profiles of type 1 and 
0.5 and 0.75. 


I type 2 configurations from/ p = 0.1 (solid lines) to 0.5 (dashed lines) and 1 (dotted lines) at/ w = 0.1 


, = 0.06,0.25, 
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N w > 0.5, influence of reducing homogeneity of microwave heating 
on i ]mw either dominates or follows the influence of varying q MW 
from f p = 0.1 to 1. In contrast, homogeneity of microwave heating 
does not alter significantly with f p in thin regime (as shown in 
Fig. 13 for N w — 0.06) and tjmw follows the variations of q MW for 
N w < 0.1 as seen in Fig. 12. In between 0.1 < N w < 0.5, changes in 
both homogeneity of power absorption and intensity of q MW 
combinedly influence ijmw as reflected in its complex variation 
patterns from f p = 0.1 to 1 within that range (see Fig. 12). 

9. Effect of f w on influence of metallic coating 

Parameter f w = A e ff/Ao alters relative magnitude of microwave 
wavelength within the sample compared to that in free space, 
which in turn alters scattering of microwaves from reactor walls via 
Eqs. (5d) and (5e). Associated changes in power absorption patterns 
for type 1 and type 2 configurations are shown in Fig. 14 for 
representative reactor dimensions of N w = 0.06, 0.25, 0.5 and 0.75, 
while the corresponding changes in q MW are shown in Fig. 15 over 
the entire range of 0.05 < JV W < 10. In both the figures (Figs. 14 and 
15), f w is varied from 0.1 (solid line) to 0.5 (dashed line) and 1 
(dotted line) while f p is kept constant at 0.1. It may be noted that 
similar to f p , f w is also bound between 0 and 1 (0 < f w < 1) as 
wavelength of microwave is maximum in vacuum or free space. 
Fig. 15 also illustrates the corresponding changes of tjmw at = 10 
(similar to Fig. 12), where other parameters used in the simulations 
are same as in previous figures (i.e. t d = 1, JV R = 0.1, /? = 0.1 and 
7R=10). 

It may be noted from Fig. 14 that unlike f p , influence of f w on 
power absorption pattern does not follow a specific trend and de¬ 
pends strongly on reactor dimension as well as reactor configura¬ 
tion. Power absorption pattern remains almost unchanged from 
fw = 0.1 to 1 in case of type 1 configuration, while f w significantly 
alters power absorption patterns for type 2 configuration for all the 
cases except for N w — 0.06 corresponding to thin regime (see 
Fig. 14). On the other hand, both the configurations interestingly 
show similar variations of q MW with suppression of resonating 
peaks from f w - 0.1 to 1 as well as higher magnitudes of q MW at 
fw = 0.5 and 1 compared to those at/ w = 0.1 over the entire inter¬ 
mediate regime except near the boundary of thin regime 
(0.1 < N w < 0.5) (see Fig. 15). Accordingly, tjmw also shows higher 
magnitudes at / w = 0.5 and 1 compared to those at f w = 0.1 over 
most of the intermediate regime along with suppression of oscil¬ 
lations as f w increases from 0.1 to 1 in Fig. 15 for both type 1 and 
type 2 configurations. These transformations of tjmw interestingly 
alter the range of reactor dimension for selection of either type 1 or 
type 2 configuration in order to enhance tjmw- At f w = 0.5 and 1, 
type 1 configuration can effectively enhance q MW or tjmw over most 



lines) and type 2 (dotted lines) configurations at f p = 0.1 and 4> R = 10. 


of the intermediate regime. This is in contrast to the case of/ w = 0.1, 
where type 1 configuration preferentially enhances q MW and hence 
tjmw only around the resonating peaks at N w = 0.5-0.25n as may be 
observed from Fig. 15. The region of preferential enhancement of 
Vmw by type 1 configuration expands to N w > 0.6 and N w > 0.45 at 
fw — 0.5 and 1, respectively. On the other hand, type 1 configuration 
no longer remains effective to enhance tjmw near N w — 0.25 at 
f w - 0.5 and 1 as may be observed from Fig. 15. Instead, type 2 
configuration becomes efficient choice to enhance tjmw till N w < 0.6 
and N w < 0.45 at/ w = 0.5 and 1 respectively in spite of lower q MW 
than type 1 configuration around N w — 0.25. This is due to the 
transformation of power absorption pattern with/ w as illustrated in 
Fig. 14, which shows that highly one sided power absorption of 
N w = 0.25 at/ w = 0.1 transforms to volumetric power absorption at 
fw — 0.5 and 1 for type 2 configuration. On the other hand, power 
absorption for type 1 configuration remains to be one sided for all 
the/ w ’s. Correspondingly, influence of homogenized heating under 
type 2 configuration at/ w = 0.5 and 1 dominate and tjmw remains to 
be higher in type 2 configuration around N w = 0.25 for/ w = 0.5 and 
1 as observed in Fig. 15. 



Fig. 14. Variations of power absorption profiles of type 
0.5 and 0.75. 


I type 2 configurations from/ w = 0.1 (solic 


s) to 0.5 (dashed lines) and 1 (dotted lines) at f p = 0.1: 
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10. Summary and conclusive remarks 

This work presents a theoretical basis for probable energy sav¬ 
ings under microwave heating for conducting an endothermic gas 
phase reaction and shows the role of metallic coating in enhancing 
the savings of energy in places where savings are not significant 
otherwise. The reaction is assumed to occur in presence of sus¬ 
ceptors within a semiinfinite batch reactor and total power con¬ 
sumption under microwave heating are compared with those 
under equivalent conventional heating. The analysis presented in 
this work are based on single mode microwave radiations, which 
can be easily achieved within any custom made waveguide. The 
entire assembly of reactor with packed susceptors can be placed in 
the test section of the waveguide, which can be further customized 
by applying a metal coated film on the opposite side of wave 
propagation to investigate the reactor dynamics with and without 
the metal coated film. It may be noted that single mode microwaves 
not only reduce complexity of the mathematical model without 
destroying the essential physics but also cause much higher power 
absorption than multi mode microwaves. The analysis has been 
performed by simulating reactor progress from detailed mass and 
energy balance equations in conjunction with Helmholtz equation 
of wave propagation by using weak formulation of Galerkin finite 
element method and Crank-Nicholson time integration scheme. In 
all the simulations, same heating rates are maintained in micro- 
wave and conventional heatings to make them comparable and 
equivalent to each other. 

This work considers probable uniform, oscillatory and exponen¬ 
tially attenuated microwave heating patterns of thin (N w < 0.1), 
intermediate [0.1 < N w < 3/(2tt/ p )] and thick (27t/V w / p or N p > 3) 
regimes and simulates them by varying reactor dimension (2L/A e tr) 
over the range 0.05 < N w < 10 for both the cases of with and 
without metallic coating. The analysis shows that microwave 
heating homogenizes reactor dynamics and can save as much as 
60% energy consumption based on regime of operation (reactor 
dimension, dielectric properties of the susceptors and transport 
limitations within the reactor). This analysis shows efficient use of 
metallic coating in enhancing savings of energy consumption under 
microwave heating for reactor dimensions around N w — 0.5n—0.25, 
n = 1,2..., where otherwise microwave heating leads to very slow 
reaction and less energy savings. For example, savings of energy 
doubles from 20% to 40% in presence of metallic coating at 
N w = 0.25 for significant mass transfer limitations within the 
reactor corresponding to = 10. Similarly, metallic coating en¬ 
hances saving of energy from 40% to 60% and 50% at N w — 0.75 and 
1.25, respectively. On the other hand, metallic coatings are shown 
to be avoided at resonating reactor dimensions without metallic 
coating corresponding to N w = 0.5n, n - 1,2. .. as well as for reactors 
in thin regime. 

This work quantifies probable savings of energy under micro- 
wave heating (i|mw) as reactor dimension varies from thin to thick 
regimes and identifies the ranges of reactors where metallic coating 
can efficiently reduce energy consumption. These results can be 
used as basis for appropriate selection of reactor configuration 
(with or without metallic coating) for a given reactor dimension as 
well as to select reactor dimension in order to maximize utilization 
of microwave energy. This work also quantifies changes in energy 
savings due to changing dielectric properties of the reacting me¬ 
dium and shows how tuning of susceptor material in conjunction 
with appropriate selection of reactor configuration (with or 
without metallic coating) can play an important role in enhancing 
energy utilization via microwave heating. 

Reduction of energy consumption is one of the important factors 
in determining sustainability of microwave heating compared to 
conventional heating to conduct chemical reaction. In that context, 


Table 1 

Literature data regarding relative savings of energy under microwave heating 
compared to conventional heating at similar operating conditions (reactor size, re¬ 
action temperature etc. 


Reference Reaction 


Heterogenous Suzuki 
Friedel-Crafts acetylatior 
Biodiesel production 


Preliminary calculations 
(Depending on 


S N Ar reaction at 140 °C 66% MARS M 

Alkylation reaction at 150 °C 59% 

Heck reaction at 140 °C 54% 

Diels-Alder reaction 11% & 71% MARS ui 


[43] FFA removal by esterification 37% 

[42] Hantzsch dihydropyridine 65% Multimode MW unit 

synthesis 


it becomes extremely essential to analyze relative energy con¬ 
sumption under microwave heating compared to that required for 
conventional heating to achieve desired yield. However, majority of 
the previous work only concentrated on the chemistry aspect of the 
process and reported the extent of enhancement of reaction for a 
variety of reacting systems under different reaction conditions. Till 
date, only a very few studies exist, which report comparative 
measurement of energy consumption under microwave heating 
and equivalent conventional heating [42-44,48,55,57], The obser¬ 
vations of these works are summarized in Table 1. It is interesting to 
note that current simulations also predict similar range of energy 
savings as measured/estimated in earlier work (especially those 
reported in the last five entries of Table 1 [42—44,48,55]) even 
though the reacting systems differ completely from the present 
setup. It is important to note that Table 1 summarizes only those 
results corresponding to similar operating parameters for micro- 
wave and conventional heating (e.g. reaction volume, reaction 
temp etc.). As noted in the previous works [42-44,48,55,57], sav¬ 
ings of energy can vary widely depending on operating parameters 
and even with the specific microwave unit used for the reactions. 
Also, reaction volume, reacting medium, mode of operation 
(continuous or batch) can significantly affect the energy savings 
among many other factors. Thus, it has been recommended by 
previous researchers to do a case-to-case basis energy consumption 
analysis to study the viability of the process. In that context, the 
present simulations may act as basic guidelines to design appro¬ 
priate reactor setup with or without metallic coating in order to 
maximize energy utilization. 
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Appendix A. Electromagnetic equations for evaluation of 
Qmw 

Propagation of microwave within a semiinfinite domain with 
dielectric packings shown in Fig. 1 can be characterized by the 
following Helmholtz equation [65] for both type 1 and type 2 
configurations: 
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and associated power absorption can be expressed a 


4ff D Peff 

In the above equations, E m is microwave induced electric field, c is 
velocity of light, e 0 is free space permittivity, /.q = c/f is wavelength 
of microwave in free space and A eff and D Pdr are effective wave¬ 
length and penetration depth of microwave within the packed 
domain of Fig. 1, respectively given by 


- y[ 2 W*& 


- 1/2 


D '--l[ 2 (77i 


where K' eff and /c" ff are effective dielectric constant and effective 
dielectric loss of the packed domain, respectively and /is frequency 
of incident radiation. Electromagnetic wave undergoes scattering at 
the interface of free space and power absorbing medium in absence 
of metallic surface/coating, which can be mathematically repre¬ 
sented by radiation boundary condition for E m [76]. Radiation 
boundary condition at the left wall is given 


The description of electromagnetic wave propagation through 
packed dielectrics under consideration is complete with appro¬ 
priate correlations K' eff and x" ff . Most widely accepted correlations 
were first developed by Fricke [77], which considers the effect of 
porosity (0) and dispersity of packings in addition to the dielectric 
properties of packing material (k' p and k") 

, , { „ ( K P + k p)f 1+ ^ + a ( K P + iK p) 0-*)] 

K eff + tK eff - -;—7—+7-T7Y7-77-> 

1-0 + [k' p + ik" J (a + 0) 

where a = 2 for spherical and 1 for cylindrical packings. 


Appendix B. Analytical expressions for absorbed power and 
their inner products 

Solving Helmholtz equation (Eq. (5c)) and associated radiation 
boundary conditions (Eqs. (5d) and (5e)), analytical expressions for 
absorbed power (<7 mw) can be written as 

16nN m fpf w [cosh 2N P (1 - z) - cos 4 ttN w (1 - z)\ 

<?mw — Ci CQS 4^^ + C2 cos h 2N P + c 3 sin 4izN w + c 4 sinh 2Np 

(B.1) 

for type 1 reactor with metallic coating and 


1 fiir/v f f f c 2 cosh 2N p( 1 - z ) + c 4 sinh 2N P (1 - z) 1 
wjpjw [ +Cl cos 4 TC ^ w (l - z) + c 3 sin 4nN m (1 - z) J 
(dj - cf)cos 4nN w + (cl + c|)cosh 2 N p - 2 CjC 3 sin 4nN w + 2 c 2 c 4 sinh 2 N p 


for both type 1 and type 2 configurations, where E 0 = ^/2/ 0 /ce 0 is 
the incident electric field. Radiation boundary condition for the 
right wall of type 2 reactor with no metallic coating is given by 


@Z = L ■■ -gjp - ij^ E m = 0 for type 2. (A.5a) 

On the other hand, metallic coating completely reflects electro¬ 
magnetic wave in type 1 reactor and associated boundary condition 
can be written as 


for type 2 reactor configuration without metallic coating, where 


coefficients ci—c 4 are given by 


Cl = 1 +/ p 2 ~fl 

(B.3a) 

c 2 = 1 +f p +/w, 

(B.3b) 

c 3 - -2/p/w, 

(B.3c) 

c 4 = 2 f w . 

(B.3d) 


Using the above expressions for absorbed power, the expression 
of (<7 mW; w„) for n + 0 can be written as 


(<?MW, Wn) 


4\/2/v, 


Cl cos 



f 2~N w / p cos 4tzN w n = 4N w ' 

sinh 2 N p j 

/pSin 47tN w 

^ 2 n+4N w 

1 + 4 N$ | 

[ 1 “ 16tc 2 N2 


4tc N w + c 2 cosh 2 N p + c 3 sin 4 tcN w + c 4 sinh 2N P 


(B.4a) 


@Z = I: E m = 0 for type 1. 


(A.5b) 
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for type 1 reactor with metallic coating and 


4V2f u 


' ( C 1 - c f) ( 


c 2 sinh 2N P + c 4 (cosh 2N P - (-1)") 

1+ W 

|" 2%N w f p (C\ cos 4nN w + C3 sin 4tzN w ) 
/ p (d sin 4 tcN w + c 3 [(-l) n - cos 47TiV w ]) 


n = 4N W 
n=t=4N m 


16tc 2 JV2 

; 47tJV w + (Cj + c|)cosh 2N p - 2c,c 3 sin 4 nN w + 2c 2 c 4 sinh 2N p 


(B.4b) 


for type 2 reactor configuration without metallic coating. Similarly, 
the expressions for (Qmw,wo) = {Qmw, 1) = q M w can be written as 


<Jmw 


4f m [sinh 2N p - / p sin 47tJV w j 

Ci cos 47 tN w + c 2 cosh 2N p + c 3 sin 4nN w + c 4 sinh 2 N p 


(B.5a) 




4 f w [c 2 sinh 2N P + c 4 (cosh 2 N p - 1) +/ p (Ci sin 4%N W + c 3 (l - cos 4ir/V w ))J 
Qmw |- c 2 _ c 2j cos 4tiW w _|_ ( c 2 + c^) cosh 2N P - 2ciC 3 sin 47tiV w + 2c 2 c 4 sinh 2N P 


(B.5b) 


for type 1 and type 2 reactor configurations, respectively. 
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